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We argue that the spontaneous breakdown of symmetry in the chirally symmetric Nambu-Jona- 
Lasinio model which was supposed to illustrate the origin of the low mass of pions in hadron physics 
does not occur due to strong fluctuations in the a - n field space. Although quarks acquire a 
constituent mass, a and n turn out to have equal heavy masses of the order of the constituent quark 
mass. 



- - - . I. INTRODUCTION 

o 

The chirally symmetric Nambu-Jona-Lasinio model |L| was the first theoretical laboratory to illustrate how light 
' pions arise from a spontaneous breakdown of chiral symmetry in hadron physics. The first realistic formulation of 
. the model which included flavored quarks, possessed chiral symmetry SU(3) x SU(3), and a spectrum of a,%,p,Ai 
mesons and their SU(3) partners, was formulated and investigated in 1976 by one of the authors ||, and has been 
[JLh ! the source of inspiration for many papers in nuclear physics in the past twenty years By eliminating the Fermi 

■ fields in favor of a pair of collective scalar and pseudoscalar fields a and 7t, as well as vector and axial vector mesons, a 
y—{ \ Ginzburg Landau like collective field action was derived. This had been studied in detail earlier as an effective action 

guaranteeing all low-energy properties of hadronic strong interactions which were known from current algebra and 
' partial conservation of the axial current (PCAC). 

In two important respects, however, the model was unsatisfactory. First it was not renormalizable in four dimensions, 

but required a momentum space cutoff A to produce finite results. Moreover, to obtain physical quantities of the 
, correct size, the cutoff had to be rather small, below one GeV, thus limiting the reliability of the predictions to very 

■ low energies. Second, if the fermions were identified with quarks, the model could not account for their confinement. 

The nonrenormalizability was removed in by replacing the four-fermion interaction by the exchange of a massive 
0^ 1 vector meson. The different attractive meson channels were obtained by a Fierz transformation. The mass of the 
0^ ' vector meson took over the role of the cutoff. The energy range of applicability was, however, not increased since the 
(-h model would still allow for free massive quarks. 

The purpose of this note is to point out a much more severe problem with the model which seems to invalidate most 
conclusions derived from it in the literature: If chiral fluctuations are taken into account in a certain nonperturbative 
approximation, the spontaneous symmetry breakdown disappears, and the zero-mass pions acquire the same mass as 
the u-mesons, both of the same order as the constituent quark mass. The nonperturbative nature of the argument 
• • ■ seems to be the reason why the phenomena has been overlooked until now. 
. £h ' Since the Nambu-Jona-Lasinio model is incapable of accounting for confinement, it gave no reason for introducing 
colored quarks. It is curious to observe that the restoration of symmetry by chiral fluctuations would offer such 
a reason, albeit with an unphysical number of colors: the physically desired spontaneous symmetry breakdown 
conclusion can only be achieved by introducing at least five identical replica of fermions. The existing three colors 
are insufficient to save the purpose of the model. 

The non-perturbative arguments used in this paper are analogous to those applied before in a discussion of the 
Gross-Neveu model [Q in 2 + e dimensions ||, where it was shown that this model has two phase transitions, one 
where quarks become massive and another one where chiral symmetry breaks spontaneously. They have also been 
applied to explain the experimental observation of two transitions in high-T c superconductors || , and to show that 
directional fluctuations in Ginzburg-Landau theories with spontaneously broken O(N) symmetry disorder the system 
before size fluctuations of the order field become relevant & . 
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II. NAMBU-JONA-LASINIO MODEL 

Let us briefly recall the relevant features of the Nambu-Jona-Lasinio model for our considerations. The model 
contains Nf quark fields ip(x), one for each flavor. Each of them may appear with N c colors, such that the total 
number of quarks is N = Nf x N c . Since the fluctuation phenomenon to be discussed will be caused by the almost 
massless modes, we may restrict ourselves to the almost massless up and down quarks. We will comment later on the 
effect of the heavier quarks. 

The Lagrangian of the model is given by || 

C = i> (0 - ™ ) V> + [OHO' + (^A a i75^) 2 ] , (1) 

where an implicit summation over a = 1, 2, 3 is assumed. A small diagonal quark mass matrix mo breaks slightly the 
SU(2) x SU (2) part of the chiral symmetry which lifts the mass of the pion to a small nonzero value. We have omitted 
the flavor symmetric vector gluon exchange used in Ref. which would have given rise, after a Fierz transformation, 
to additional vector and axial vector interactions, which would not influence the chiral fluctuations to be investigated 
here. Thus we use the original nonrenormalizable interaction corresponding to an infinite vector gluon mass. The 
coupling constant in (|l|) is defined with the number of colors N c in the denominator, to allow for a finite N c — > oo 
limit of the model at a fixed go. The 2 x 2-dimensional matrices A Q /2, with a = 1, . . . , 3 generate the fundamental 
representation of flavor SU(2), and are normalized by tr(A a A{,) = 2<5 a {,. 

Via a Hubbard- Stratonovich transformation, the Lagrangian (^) is converted into a theory of collective scalar and 
pseudoscalar fields a and 7r a . Defining the propagator in the presence of the meson fields 

G^- B . . , (2) 

and adding external quark sources 7], fj, one can integrate out the quark fields from the corresponding Lagrangian. 
Summing over colors, the generating functional of the Green functions takes the well-known form 



J VaVnexp^iNc -iTr' In iG^ 1 - J d D x {a 2 + ir 2 a ) + iN- 1 J d D xd D yf]Gr] j 



(3) 



The term inside the brackets is the collective field action A[a, tt], whereas the symbol Tr' includes both, the functional 
spacetime "index" x, and the internal trace over spin and flavor indices: Tr' = j d D xtx^tXf. 

By extremizing A[a, n] at zero sources 77,77, we obtain the field equation for the collective field (tr, n a ): 



tr 7 tr/ 



G(x, x) 



1 



90 \Ka(x) 



(4) 



For constant fields, this equation becomes a gap equation. Its solutions will be marked by a superscript "s" for 
"stationary phase approximation" . From now on, unless explicitly stated, we shall consider the model with zero mass, 
m = 0. The stationary pseudoscalar solutions 71^ can always be chosen to be vanishing, while the scalar solutions 
can be a s = 0, or er s = p a . In the first case, the ground state is chirally symmetric, in the second the symmetry is 
spontaneously broken. This is the state of physical interest whose stability will now be discussed. 

III. EFFECTIVE POTENTIAL AND GAP EQUATION 

In the limit N c — > 00, the generating functional is given exactly by the extremal field configurations, which will be 
parameterized as (o~ s (x), 7r^(x)) = (po(x), 0). The system has an effective action per quark 

Ei^-M) = -TTr'lnzG; 1 -±J d D x Po 2 + ± J dPzV.iCtfV., (5) 

where ^ = iG Po r] is the expectation value (ip) of the quark field, and G Po its propagator 

Gpo = t^—. (6) 
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This shows that the solution of the gap equation with po ^ describes quarks with a nonzero mass M = po, which 
has been generated by the spontaneous symmetry breakdown, and is referred to as the constituent quark mass. In the 
present approximation of zero bare mass mo, the constituent quark mass is about equal to 300 MeV for up and down 
quarks (see the discussion in Refs. ||||). In cither case, the Green function (||) in the stationary field is diagonal in 
flavor space. 

In the absence of external quark sources, the ground state expectation value of a fermion field is always zero, and 
the expectation value po(x) is constant, so that (||) reduces into 



r(p) 

Nr. 



-iTr'lniGp 1 



1 

Wo 



d D xp 2 , 



(J) 



where we have allowed the fields a and 7r a to be nonextremal, defining a 2 

the extremum. This is determined, after a Wick rotation to euclidean momenta p E with po — ipE,o, id v p 



tt 2 = P 2 , and reserving the notation p\ for 

D - - -d D p El 



P 



-p E , by the gap equation 



1 = 2 x 2^/ 2 f d ° PE 
9o 



(2n)Dpl + p 2 



(8) 



We have divided the two sides of the gap equation by a common factor po , since we want to study the spontaneously 
broken phase. 

The gap equation must be regularized, which may be done in many ways. Here, we shall use two methods: analytic 
continuation in the dimension D, and a cutoff A in momentum space. The former is mathematically more elegant and 
has the advantage of relating the properties in four dimensions to those in 2 + e. It has, however, some unphysical 
properties which require special attention, as we shall see. Such problems are absent in a cutoff regularization scheme, 
which exhibits clearly the physical divergences caused by the infinite number of degrees of freedom of the field system. 
Factorizing the integral in (ph into direction and size of the momentum p E , we bring the gap equation to the form 



1 

90 



2po 



D-2 



r(i - d/2) 

(2n) D / 2 ' 



(9) 



Denoting by f2 the £>-dimensional volume J d D x, the volume density v(p) = —T(p)/Q of the effective action (0) is the 
effective potential per quark. Performing the internal traces, and subtracting a divergent constant term associated 
with the chirally symmetric state with po = 0, we obtain the condensation energy in euclidean space: 



Av(p) 



Nc 
2 



9o 



D (2ir) D / 2 



r(i - d/2) 



(10) 



In an even number of dimensions D, both the gap equation (|J) and the effective potential (|h]) are divergent, due to 
a pole in the factor T(l — D/2). Introducing the diverging parameter b € — 2Y{1 — D/2)/{D(2tt) d / 2 ], we can rewrite 
the gap equation and effective potential in the more compact form as 



1 

9o 



Av(p) = ^ 



Dpo u -'b. 
' 1 



-p 2 - 2p D b t 
9o 



(11) 
(12) 



In the more physical regularization with a cutoff A in momentum space, these expressions look more complicated: 

1 2 

9o 



(2ny 



A 



A 2 

po z In ( 1 + — 

Po 2 



. . . N c f 1 2 
2 190 



(2ny 



2A2 



p 2 A 



A 4 / p 2 \ p 4 , / A 
— In I 1 + ] - f^r In I 1 + — 



(13) 
(14) 



The results (|llj) and (|12|) of the analytic regularization scheme can be mapped roughly into the cutoff results 
( |l3| ) and ( |l4| ) if we recall the special property of dimensional regularization that all integrals over pure momentum 
powers vanish identically: / d D k (k) a = (Velt man's rule). Thus, arbitrary pure powers of the cutoff A a+D have 
no counterpart in dimensional regularization. Only logarithmic divergences can be related to diverging pole terms 
\/e — > oo for e — > 0. It is therefore inconsistent to relate e to A by setting T(e/2 — 1) « A 2 /pg, as proposed by 
Krewald and Nakayama Q. Only the logarithmic divergence in ( |l3| ) can be mapped to the small-e divergence in 
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(pi]), setting r(e/2 — 1) « — ln(l + A 2 //} 2 ,). With their inconsistent identification, Krewald and Nakayama matched 
A by an e > 2 which lies in the wrong region D < 2, the physically relevant range being D E (2 + e, 4 — e). 
Note that the matching of the logarithm at the level of the effective potential leads to the properly matched gap 
equation, thus having circumvented the unphysical properties of the analytic regularization. The free use of this 
scheme in renormalizable field theories relics on the fact that all infinities are eventually absorbed in unobservable 
bare quantities, such that the artificial zeros of the integrals over pure powers of momenta cannot produce problems. 
In nonrenormalizable theories, on the other hand, only a cutoff (or a related Pauli-Villars regularization) is physical, 
and analytic regularization must be treated with caution. This is seen even more dramatically in integrals which do 
not have logarithmic infinities. For example the condensation energy (10) in D = 3 dimensions would be a finite 
negative number in analytic regularization, while being a linearly divergent positive function of the cutoff. 



IV. CHIRAL FLUCTUATIONS 

Since the physical number of quarks N c is finite, the fields perform fluctuations of magnitude l/y/N c around their 
extremal value. As long as iV c can be considered as a large number, the deviation from the extremal field configuration 
(a',Tr' a ) = (a — po,Tr a ) are small, and the action can be expanded in powers of (a',Tr' a ). The quadratic terms in this 
expansion define the propagators of the collective fields (a',Tr' a ). The higher expansion terms of the trace of the 
logarithm in (|^) define the interactions. With this decomposition, the inverse of the quark propagator (0) can be 
decomposed into a constant and a fluctuating part, setting iG^ 1 — iG~^ — (a' + ij5\ a Tr' a ), with G Po of Eq. (|6|). Then 
we have 

Tr In ^G- 1 = Tr In iG" 1 + Tr In [1 + iG Po (a' + i-ysXaK)} • (15) 

An expansion of the last term up to the second order in the fields gives an approximate partition function (with 
Z = exp [— Q E N c v(po)] and fl E is the euclidean volume) 



Z = Z J VaVirexp (iN c j^Tr' [iG Po (a' + i^XaK)] 2 d ' ^ ^ + ) ■ 



(16) 



The functional matrix between the fields in the exponent gives us directly the inverse of the desired collective free 
field propagators G a > , G^' ■ In momentum space, we identify 

M°'y] = \ J d D q[K( q )G- 1 n' a (-q) + a'( q )G- 1 a'^ q )], (17) 



where 



G-\ = 2 x 2^N C f 1 dy [ & +PE1E + {2po*,0) 

Jo J ^) D [{ql + 2p EqE )y+pl+p^ 



(18) 



In this expression, the gap equation (||) has been used to eliminate the term l/.go- The notation (2/9o 2 ,0) indicates 
that only the equation for a contains an extra term 2po 2 • 

In four spacetime dimensions, the integral evaluated in dimensional regularization reduces to q E /2 for the pseu- 
doscalars, and to (q E + 4po 2 )/2 for the scalars, both with a diverging coefficient. The first leads to a zero mass for 
pions as a manifestation of Goldstonc's theorem, the second to a mass equal to twice the constituent quark mass for 
the tr-mesons. For a finite result, the integrals must be regularized. In D = 4 — e dimensions, the inverse euclidean 
propagator is seen to start out for small q E like 

*N c (l- I) = Z^(p )q% + 0{q%). (19) 



with the same b e as defined above Eq. (11). If the theory is regularized with a cutoff A in momentum space, this 
becomes 



(27T) 2 



A 2 \ A 2 

In 1 



po 2 J A 2 + p 2 



q E = Z^(p )q E . (20) 



In the right-hand part of the two equations, the factors in front of q E have been identified as the wave function 
renormalization constants Z^po) of the pion field in the two regularization schemes. 



5 



As a consequence of the spontaneous symmetry breakdown, the fluctuations of the pseudoscalar fields are massless. 
These fields appear in the x-space version of the action ( |l7|) in a pure gradient form 

MA = { J d D x{[dKM 2 }, (21) 

with (3 = Z n . Due to chiral symmetry, this gradient action can be extended to the gradient action of an arbitrary 
field (a, 7r ). Introducing the directional unit vector fields n, = (a',n' a ) = (a',n' a )/p, we find: 

Ao[ni] = ^p-p 2 J d D x(d ni ) 2 , i = l,...,N n , (22) 

with N n = 4 and 

= Z n (p). (23) 

This chirally invariant action describes the massless pions with all multipion interactions]^] The prefactor (3 is called 
the stiffness of the directional fluctuations P]7],|ll)|~|l2| . In analytic regularization, the result ( |l9| ) shows that the 
stiffness of pion fluctuations in D—2 dimensions becomes 

2ttPo 

thus coinciding with the stiffness calculated in Ref. Q in the Gross-Neveu model (which contained a factor TV to be 
identified with the present Nf x N c — 2N C ). 

With the more physical cutoff regularization in D = 4 dimensions, the stiffness of directional fluctuations is 



1+1 - 

Po 



A N 2 



A 2 



Po 2 + A 2 



(25) 



This is the crucial quantity leading to our fatal conclusions for the restoration of chiral symmetry. The stiffness fl2q ) 
is far too small to let the directional field settle in a certain direction, required for spontaneous symmetry breakdown. 
The disordering effect of phase fluctuations is well-known from many model studies of the 0(4)-symmetric Heisenberg 
model on a lattice. High-temperature expansions and Monte Carlo simulations have shown that there exists a critical 
stiffness below which the system goes over into a disordered state. 

For an analytic estimate of the critical stiffness, we relax the unit vector constraint for the vectors m in ( p2| ) by 
introducing an additional field A(x) playing the role of a Lagrange multiplier. The rii-fields can then be integrated 
out in the partition function, leading to an action 

S = ^Trln [-d 2 + X(x)} - P{p 2 )p 2 J d D X ^L, (26) 

where Tr denotes the functional trace (the summation over the fields component has already been performed). For a 
large number N n of components, the fluctuations are suppressed, and the field A (a;) becomes a constant satisfying a 
second gap equation 



B=^f dDk 1 (27) 



If there is a nonzero solution A =/= 0, this will play the role of a square mass of the n^-fluctuations, and represents an 
order parameter in the directional phase transition. The model has a phase transition at a critical stiffness 



1 Only two approximations are involved: the first one consists in freezing the size p of the fluctuations. The second one neglects 
corrections due to the finiteness of the sigma mass. The latter corrections are expected to be of the order /^/(4Af 2 ) ~ 3%. 
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For a smaller stiffness, the phase fluctuations are so violent that the system goes into a disordered phase with A 7^ 
giving all fields a nonzero square mass A. Since the fields rii are the normalized a and 7r a fields of the model, this 
determines an equal nonzero square mass of a and ir a mesons, and thus a restoration of chiral symmetry. 

Note that the quarks are still massive: their constituent mass is a consequence of the formation of the pairs, which 
are strongly bound for small N c . The phase transition taking place at the critical value of the stiffness, on the other 
hand, is related to the Bose-Einstein condensation of the pairs. At small N c , the two processes are widely separated. 
This separation of the two transitions (pair formation and pair condensation) can be judged by the simple fluctuation 
criterion in Ref. 0. 

In our model, the number N n is equal to four, which is not very large. Fortunately, Monte Carlo studies of the 
model Jf^-[T^| have shown that N n — 4 is large enough to ensure the existence of the transition and the quantitative 
reliability of the theoretical estimate of the critical stiffness (p8|). From an evaluation of (28h on a lattice, and a 
comparison with Monte Carlo studies, we estimate that the critical stiffness obtained from ( 28| ) is correct to within 
less than 2% or 6% The same maximal error is expected if we work in the continuum using a 

momentum cutoff scheme. 

For N n — 4 and a cutoff in the integral over pion momenta, the critical stiffness is given by 



/3c 



4 A 2 



I671" 2 p 2 

By comparing this with the stiffness of the model in (p5[), we find 



(29) 



N c = 



K 
A 









jln 


I + 











(A/Po) 2 
1 + (A/po) 2 



(30) 



This equation determines the number N c of identical quarks which is necessary to produce a large enough stiffness 
j3 to prevent the restoration of chiral symmetry. Only if the number of colors exceeds this critical value, will the 
model possess a phase in which the pion is a massless Goldstone boson, and a a meson with a mass twice as large 
as that of the constituent quarks. The critical number (30) is plotted as the solid curve in Fig. for K = A. We 
see that N c = 5 would be the smallest allowed value. This number, however, is incompatible with color SU(3). This 
suggests that the Nambu-Jona-Lasinio model always remains in the symmetric phase, due to chiral fluctuations. It 
can therefore not be used to describe the chiral symmetry breakdown of quark physics, as has been claimed by many 
publications, which have appeared in particular in nuclear physics ||. 

Can this conclusion be avoided by a different choice of parameters? To obtain a critical value smaller than N c = 3 
would require a pionic cutoff A^ < 0.8A. However, the cutoff cannot be chosen at will. Let us study the cutoff 
dependence more precisely. For this, we refine the previous crude estimate d28|), ( p9| ) of the critical stiffness, which 
will henceforth be called Approximation 1, by taking better account of the shorter wavelength fluctuations, replacing 
the action (E2T) by 



Ai[m} 



P_ 

2 



d D xn l {x)G- 1 (~d 2 )n i (x), 



(31) 



with G^T 1 (— d 2 ) from Eq. ([l8|). This exchanges 1/fc 2 in Eq. (28) by the full pion propagator G w (fc 2 )/zi A ^ associated 
with the action (|3l|). The cutoff A T makes the integral over pion momenta finite. Its size is fixed by physical 
considerations. The pion fields in the symmetry-broken phase are composite, and will certainly not be defined over 
length scales much shorter than the inverse binding energy of the pair wave function, which is equal to 2M = 2po. 
Thus we perform the integral in the modified equation ( p8| ) up to the cutoff 4Af 2 . This is Approximation 2, yielding 
the solid curve in Fig. ||. 

The phase with broken symmetry for three colors would be reached only if the quark loop integration is cut off at 
A 2 > 11M 2 . Such a large value, however, is incompatible with the experimental value of the pion decay constant 
fn ~ 0.093 which is given, in the large- N c limit of the model, by f^/M = Z X I 2 {M). For typical estimates of 
constituent quark masses M e (300,400) MeV @, we find that A 2 /M 2 should lie in the range (3.3,7.3), the highest 
value corresponding to the lowest possible mass 300 MeV. 

The above study has given us only the critical point, where the pion mass goes to zero. We can do more and 
determine the common nonzero square masses to 2 = m 2 = A of a and 7r a -fields in the phase of restored chiral 
symmetry. This is the subject of the next section. 
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V. MESON MASSES 



The chiral fluctuations give rise to a change of the effective potential. They add to Av(p) in Eq. (|l4j) an additional 
energy coming from the stationary point of the action (Eq) at a constant A(x) — A: 



XZ oP 2 + 



16tt 2 



m[<z| + A], 



A' 2 v(p,X) = --XZ(p)p 2 + ^ 



, : N n /• u dq 2 E q 2 E u ^, , 



16tt 2 



ln[G- i (4)/Z(p) + A] ! 



(32) 
(33) 



for Apprs. 1 and 2, respectively, where the latter has — d 2 replaced by G 7r 1 (— d 2 )/Z{p). Extremizing Av(p) + 
A'^ 2 ) v (Pi X) yields two coupled gap equations replacing the independent gap equations ( |l3| ) and (p7f). Introduc- 
ing the reduced quantities Z(x) = In (l + a; -1 ) — (1 + x) -1 , and x = p 2 /A 2 , y = A/A 2 , we have for Appr. 1: 



x In (l + x r ) + ^— [sZ(i)] =xln(l + a; 



( (x ) -■" ,n 



(34) 
(35) 



For Appr. 2, the coupled gap equations are more complicated since the full g 2 -dependence of Z n has to be taken into 
account. They read 



r-(A„/A) 2 



Z(x ) 



^(l + x- 1 )-^ 



Z(k 2 ,x a ) 

(A. /A) 2 



d 

dxQ 



' Z(k 2 ,x ) 
Z(x ) 

k 4 dk 2 



yd_ 

2 dx 



[xZ(x)] 



N c xZ(x) 



SN C |7 k 2 [Z(k 2 ,x)/Z(x)] +ydx 
X,,} HAtt/a ) 2 k 2 dk 2 

o k 2 [Z(k 2 ,x)/Z(x)]+y' 



Z(k 2 ,x) 



Z(x) 



where Z(k 2 ,x) is taken from the pion propagator (pq), and is given by 



Z(k 2 ,x) 



Nc_ 
2tt 2 



p 2 dp 2 



(1 — z)dz 



[p 2 + k 2 z(l - z) + '. 



(36) 
(37) 

(38) 



There is no need to write down the lengthy analytic solution to this double integral. The coupling constant go has 
been eliminated in favor of the reduced mass xq which characterizes the model uniquely above , where y = and 
thus A = 0. In that case, Eq. (34) of Appr. 1 reduces to (^3|). Equations ( j35"l) and (p7|), on the other hand, determine 
the common square mass A of a and 7r a as a function of N c , which begins developing for N c < N". Note that going 
from Appr. 2 to 1 corresponds to using a momentum-independent pion normalization Z{k 2 ,x) = Z(x). This makes 
© and @ coincide with @) and ©. 

The solutions of (^4|)-(|3^) are plotted in Figs. Figs. are for Appr. 1, restricted to the case = A. 

Qualitatively, the pictures remain the same for different ratios A^/A. Quantitatively, there is only a shift in the 
critical number of color (solid curve of Fig. [TJ) to A^ r = 3 as A^/A is lowered to 0.8, while it increases above the 
given curve if A^/A > 1. This is due to the fact that at the critical point corresponding to A = 0, one sees from 
Eq. © (or from Eq. (||) with y = 0) that A^ cr oc (A w /A) 2 . The dashed curves of Figs. |l| and || are explained in the 
corresponding legends. Here we only remark that the shape of the dashed curves in Fig. [l] can be understood from 
the gap equations (34) and ( |35| ) without solving them, because xZ(x) is maximal at the minimum of N". 

Figures and ^ correspond to Appr. 2, in which the full momentum dependence for the pion normalization constant 
is taken into account, and in which the pionic cutoff is A 2 = 4M 2 , for which we get the ratio (A^/A) 2 = 4x . The 
solid curve in Fig. |^ gives the critical number of color in this particular case. Although the conclusion is not as strong 
as in Appr. 1, our result concerning the lack of breaking of chiral symmetry is robust, since the crossing with the 
line N c = 3 takes place at a cutoff A 2 /M 2 > 11, which lies outside of the admissible range (3.3,7.3) implied by the 
physical value of the pion decay constant = 93 MeV, as discussed at the end of the previous section. 
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Finally, we give in Fig. || the stiffness as a function of the number of color for Appr. 1. The three curves depend so 
weakly on p that they seem to coincide. To make the /?o~dependence visible, we have plotted an extra dotted curve 
for a very small value po — 0.224A (dotted). 

Let us emphasize that these conclusions cannot be reached in the dimensional regularization scheme since, as 
explained at the end of Section 3, the integral in (psj ) determining the critical stiffness vanishes. Here the unphysical 
nature of dimensional regularization makes its application impossible. 

Before concluding, let us also remark that the cutoff chosen in Appr. 2 is completely different from that in Appr. 1, 
where the ratio of cutoffs is a constant. In Appr. 2, the ratio of cutoffs is a function of xq: A 2 /A 2 = 4xq. If we had 
taken the cutoff in the same way as in Appr. 1 (A 2 /A 2 = 1), the curve giving the critical number of colors would also 
have had the same shape as in Appr. 1, although the integration would have been much more involved: the minimum 
number of color would then be 5.2, whatever the value of A 2 /pq, a value which is even higher than in Appr. 1. We 
see that Appr. 2 as presented above, with the physically motivated cutoff A 2 = 4M 2 , gives then the lowest critical 
number of colors. 

Our conclusions were derived from a study of only the a, 7r fields. The inclusion of other flavors does not help 
preventing the restoration since the associated pseudoscalar mesons are quite massive, making their fluctuations 
irrelevant to the described phenomenon. 



VI. CONCLUSION 



We have shown that within a certain nonperturbative approximation, the Nambu-Jona-Lasinio model does not 
really display the spontaneous symmetry breakdown for whose illustration it was constructed. The fluctuations 
of a- and 7r a -fields restore chiral symmetry and make a and tt equally massive. If our conclusion survives more 
refined approximations, this would invalidate a large number of publications, especially in nuclear physics, which 
have been based on the existence of a symmetry-broken ground state of the model. In particular, all studies of the 
temperature dependence of the symmetry-broken state would deal with nonexisting objects, thus calling for further 
investigations. Finally, we note that our no-go result for the Nambu-Jona-Lasinio model does not imply problems 
with the effective-action approach to chiral dynamics. Certainly, there exists an effective chiral action for the meson 
sector of quantum chromodynamics which does contain almost massless pions for N c = 3. It is only the Nambu- 
Jona-Lasinio model as it stands which is incapable of describing these for such a low number of colors. In fact, a 
recent paper fll7t prompted by a first version of our preprint points out that an extension of the Nambu-Jona-Lasinio 
model by interactions involving higher-dimensional operators is not subject to our no-go theorem. Another escape is 
possible by adding gradient and quartic interaction terms for a- and 7r-fields to the initial action, thus extending the 
Nambu-Jona-Lasinio model to a linear sigma model llq] . 
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FIG. 1. Solid curve shows Approximation 1 to critical num- 
ber of colors TV" as a function of the extremal value of 



p = Vo" 2 + 7i"a, above which chiral symmetry is restored. The 
dashed curves indicate the solutions to the two gap equa- 
tions ( pS4| ) and ( pS5| ) for three different values of the constituent 
quark mass po = M in the symmetry-broken phase above N£ T ■ 
The three quark masses lie above (po > p*), below (po < p*), 
and at p* = M* « V0.46A where A 7 "" takes the minimal 
value 4.62, with a constituent quark mass above N" of 0.678A 
(short-dashed curve). The medium-dashed curve corresponds 
to a constituent quark mass 0.479A, and the long-dashed to 
1.342A. 
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FIG. 3. Same plot as in Fig. (|fj), but for Approximation 2, 
the dashed curves indicating the solutions of the two gap equa- 
tions © and ©. The three mass values are now po > p* , 
po < p*, po = p* , with p* = M* « V0.092A, corresonding to 
the minimal critical number of colors iV" = 3, implying a con- 
stituent quark mass above JV" = 3 of 0.303A (short-dashed 
curve) . The medium-dashed curve are for a constituent quark 
mass 0.447A, and the long-dashed for 0.224A. 
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FIG. 2. Common square masses m 2 = m 2 = A as a func- 
tion of N c in Approximation 1. The three curves start at 
different critical values A^ 1 which can be read off Fig. [j]. 
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FIG. 4. Common square masses m 2 = m 2 = A as a func- 
tion of N c in Approximation 2. The three curves start at 
different critical values N^ T which can be read off Fig. H. 
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FIG. 5. Approximation 1. Reduced stiffness as a function 
of N c /N" . The three curves below N™ cannot be distin- 
guished on this scale. The dotted curve corresponds to an 
extra low value of po, just to show that below N^ T the curves 
deviate from a straight line. 



